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A nonlinear Fokker-Planck equation is obtained in the continuous limit of a one-dimensional lattice 
with an energy landscape of wells and barriers. Interaction is possible among particles in the same 
energy well. A parameter 7 , related to the barrier’s heights, is introduced. Its value is determinant 
for the functional dependence of the mobility and diffusion coefficient on particle concentration, but 
has no influence on the equilibrium solution. A relation between the mean field potential and the 
microscopic interaction energy is derived. The results are illustrated with classical particles with 
interactions that reproduce fermion and boson statistics. 

PACS numbers: 05.40.-a, 05.60.-k, 66.10.Cg, 05.10.Gg 


I. INTRODUCTION 

The complete description of a large number of interact¬ 
ing particles (many-body system) is a problem that usu¬ 
ally exceeds numerical or analytical capabilities. Mean 
field theory allows the description of such a system by 
analyzing the behaviour of only one particle subjeted to 
an average force produced by the rest of the particles. 
The mean field force depends on statistical properties 
of the rest of the particles that, for self consistency, are 
equal to the properties of the particle considered before 
[TJ p. 131]. This is the origin of the nonlinear char¬ 
acter of the resulting description [Up. 3]. The non- 
linearities reflect the interactions between particles. For 
example, Kaniadakis and Quarati |3] introduced a non¬ 
linear Fokker-Planck equation for the particle density in 
a classical system that reproduce quantum statistics. An 
appropriate choice of the transition probabilities, that 
depend on the particle density, gives rise to Fermi-Dirac 
or Bose-Einstein distributions in equilibrium. Quantum 
effects can be reproduced in a classical system with an 
interaction potential, also called statistical potential jH 
p. 124]. For example, the Pauli exclusion principle for 
fermions is analogous to a potential that becomes infinite 
when two particles occupy the same state. This hard¬ 
core interaction is the cause of a nonlinear term in the 
corresponding diffusion equation, see E], H p. 280] or 
Eq. (10) in |S]. A general approach for the derivation 
of a nonlinear Fokker-Planck equation starting form the 
Master equation can be found in EE!- More examples 
can be found in E]. 

Many researchers have devoted attention to the prob¬ 
lem of diffusion with interaction, in many cases motivated 
by the seminal work of Batchelor [8] . The following list is 
a limited and partial sample of related references: mia. 

The problem that we wish to address, in a quite general 
perspective, is collective diffusion of interacting particles. 
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We use a one dimensional lattice, but the results can be 
easily extended to higher dimensions. In our analysis, 
only the collective or transport diffusion coefficient is in¬ 
volved (we do not analyze here single particle diffusion, 
or self-diffusion of tagged particles). 

The restriction on the interaction is that it is local: 
it takes effect only among particles in the same site. In 
other words, the interaction range is smaller than the lat¬ 
tice spacing. We apply a mean field approximation as¬ 
suming that the evolution of the system can be obtained 
analyzing the behavior of only one particle. We consider 
that this particle is subjected to a mean field potential 
Vi that depends on the number of particles in the site, 
rii, where i is the lattice site index. The variation of the 
particle number is smooth so that the continuous limit 
can be applied, and a nonlinear Fokker-Planck equation 
is obtained. The equilibrium solution is completely de¬ 
termined by the mean field interaction potential and, if 
present, an external potential. This is not the case for 
the non-equilibrium behavior. As we show in the next 
sections, it depends on an additional parameter, 7, that 
determines if the transition probability between neigh¬ 
boring sites depends on the potential in the source site, 
in the target site, or on a mixture of both potentials. 

II. TRANSITION PROBABILITIES 

In the one-particle picture, the energy associated to a 
particle in site i is given by 

E, = V + C/„ (1) 

where Vi, the internal or mean field potential, is an ab¬ 
breviation of Vijii), and Ui is an external potential. We 
can consider Vi as the function of rij that, in equilibrium, 
satisfies the relation rieq^ oc exp{—/^[^(neqy) + Ui]}, i.e., 
it is a one-particle effective potential that satisfies Boltz¬ 
mann statistics in equilibrium. 

The interaction energy of a configuration of rii particles 
is given by a function $(ni). The relevant problem of 
determining the relation between the mean field potential 
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Vi and the interaction energy $(ni) is addressed in Sect. 

, The interaction is local in the sense that particles 
interact only when they are in the same site. 

In the one-particle picture, the detailed balance condi¬ 
tion gives a relation between transition probabilities: 

= ( 2 ) 

where Wiy+i is the transition probability from site i to 
site j-l-1. Using the words of Derrida [24], it is a “straight¬ 
forward generalization” to consider that (§ also holds for 
the transition probabilities out of equilibrium (it is clear 
that this does not mean that detailed balance holds out 
of equilibrium, since in general u, ^ const x 

Detailed balance is not enough to determine the tran¬ 
sition probabilities. In order to obtain an expression for 
them, we assume that they can be written as a combi¬ 
nation of exponentials of Vi, U+i, Ui and Ui^i. In the 
paragraphs below we present a physical interpretation of 
the result that also serves as a justification of this as¬ 
sumption. We have 

Wi^i+i = P exp [-/3 (7 Vi+i +^'Vi + a Ui+i + a Ui )], 

( 3 ) 

where 7 , 7 ', a, a' and P may depend on position i; 
the subindex is omitted in order to lighten the notation. 
We assume that the system has an inversion symmetry. 
It implies that the reversed transition probability is ob¬ 
tained from ([^ by exchanging i i -\-l: 

Wi+i^i = P exp [-/3( 7 Vi + 7 ' Vi+i +aUi+a’ Ui+i )], 

( 4 ) 

It can be shown that the combination that fulfills detailed 
balance is 7 ' = 7 — 1 and a' = a — 1. Then, 

Wii+i = 

A physical interpretation of 7 and a parameters is re¬ 
lated to the barriers height between neighboring sites. Le 
us consider that particles can occupy discrete sites in a 
lattice with a potential that has a continuous shape, as 
shown in Fig. Between sites i and i + 1, the potential 
has a maximum value of i?i+i/ 2 . A particle that jumps 
from z to J -|-1 has to overcome a barrier Bi_^.i /2 — Ei, and 
the transition probability is given by 

IU,,,+i ( 6 ) 

where v is the number of jump attempts per unit time. 
The potential maxima depend, on one hand, on the char¬ 
acteristics of the substratum, that contributes with a 
term C. They also may depend on the values of the 
interaction potential at both sides of the barrier. This 
dependence is symmetrical (the same for both sides) and 
its influence is represented by a term 7 (U -I- U+i), the 
choice of the name of 7 parameter advances that it is the 
same than the one introduced in ([^. Finally, the external 
potential also has an influence on the barriers. It is given 
by an additional parameter rj G ( 0 , 1 ): TfUi -I- (1 — r/)Ui+i. 
The influence of the external potential on the barriers is 


7 = 0 



7 = 1/2 




FIG. 1. Energy landscape of wells (Pi) an barriers (-Bi+i/ 2 )- 
From to to bottom: 7 = 0, 1/2 and 1 Vi. In this illustration, 
C is constant and Ui — 0 Vi. 

always present. Parameter rj is an interpolation factor 
between sites i and z -|- 1 for this influence. Taking these 
arguments together, the potential maximum is 

Bi+l/2 = C' + 7(U + U+i)+ryC/, + (l-z7)C/,+i. (7) 

Replacing in and considering that P = ve~^^, we 
obtain 

Wii+i = . (8) 

Let us compare this with Eq. (§. Since the values of 
Ui and are, in principle, arbitrary, we have that 

a— 1 = —l-l-zy and a = 1 — 77 and, therefore, a = rj = 
1 /2. This result is actually a consequence of the inversion 
symmetry assumed in Finally, the expression for the 
transition probability is 

Wi^i+i = ( 9 ) 

with AU = C/j+i — Ui. 
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The previous arguments allow us to construct an en¬ 
ergy landscape that gives physical meaning to the param¬ 
eters and variables involved in the transition probability. 
In this picture we are using the Hanggi [5^ interpretation 
of Lagevin equation with multiplicative noise, or barrier 
model [26], for which the particle current, without exter¬ 
nal potential and in a continuous space, is J = —D{x)^ 
(for other interpretations, the spatial dependent diffusion 
coefficient is inside the space derivative | 2 S|). 

In principle, 7 can take any real value. We conjecture 
that, from a physical point of view, the pertinent values 
of 7 are in the range [0,1]. Typical situations are shown 
in Fig. For 7 = 0, the transition probability depends 
on the origin potential: 

= Pe-0(-y^+^u/2) ^ 

For 7 = 1 , the transition depends on the target potential: 


does not depend on 7 . In the absence of interaction, we 
recover the Einstein relation: D = . 

The resulting Fokker-Planck equation 


dn d f dU ^ dn\ 


(18) 


is nonlinear because of the dependence of ^ and D on 
n and V(n). It is a free energy Fokker-Planck equation 
with Boltzmann statistics, see chapter 5 in |5]. If we 
consider the notation used in Eq. (5.4) of reference |2|, 
we find the following correspondence: n —> P, PU —)■ 
Uq/Q, PV —>■ and Dq —>• Q. The main novelty of 

our approach is the introduction of 7 parameter —that 
plays a relevant role in the system’s dynamics— and its 
physical interpretation in the energy landscape. 

Identifying the zero current state with equilibrium, it 
is easy to see that the equilibrium concentration is rieq oc 
Q-P(v+u)^ The proportionality constant can be written 
as 


And for 7 = 1/2 we have an intermediate case that is 
a frequent choice in Monte Carlo (MC) simulations of 
diffusion processes; the transition depends on the energy 
difference between target and origin potentials: 

(y = l/ 2 ), ( 12 ) 

with AV = Vi+i — Vi- 


= (19) 

where we can identify fic with the chemical potential. As 
expected, 7 plays no role in the equilibrium solution. Its 
influence is present in the dynamics through the depen¬ 
dence of ^ and P on 7 . Let us also note that the 7 pa¬ 
rameter has influence only when an interaction potential 
is present. 


III. CURRENT AND FOKKER-PLANCK 
EQUATION 


The current J between sites i and f + 1 is 




(13) 


We replace in ( |13[ ). Now we turn to a continu¬ 
ous description in which rii is replaced by n{x), Vi by 
V(n{x)), and Ui by U{x), with x = ai; we approxi¬ 
mate {rii+i - n^)/a ~ {Vt+i - Vi)/a ^ and 

{Ui+i — Ui)/a ~ We call Dq = PaP the free diffusion 
coefficient. After some algebra, we obtain 


, dU ^ dn 

Ja = -n — n- D —, 
ox ox 


(14) 


where the mobility fj, and the diffusion coefficient D are 
= PDoe-^P''^-'^'^^ (15) 

D = Dq 


In the previous derivation the validity of the Ginzburg 
criterion for a mean field theory is assumed: the fluctu¬ 
ations are small enough so that ((n/) — {ni)^)/{n'f) 1. 

The relation between mobility and diffusion coefficient 


IV. THE MEAN FIELD POTENTIAL 


Instead of speaking about the energy of one particle, 
as in 0 ^ let us consider the energy of the conhguration 
of Ui particles: 


Ei = ^{ui) + mUi, ( 20 ) 

where ^{ui) is the interaction energy of the configuration 
of Hi particles. The question that lies at the core of a 
mean field theory is: what is the relation between V and 
$? 

In order to answer this question we have to appeal to 
the microscopic description given by the grand partition 
function for classical particles. Since there is not inter¬ 
action between different lattice sites, it can be written 
as 


E=nz. 


i 


( 21 ) 


with 


= 



ml 



^ —Sini—0^(ni) 

nf 


( 22 ) 


where, for simplicity, we introduced 5i = P{Ui — ^c)- We 
can obtain the equilibrium mean value (rii) with 




(Ui) 


d\nZi 

~w~ 


D = ftP 


(23) 
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and express this result as a function of Si (for a given 
value of /?): (rii) = fp[Si). This means that we can 
consider Si as a parameter that allows us to scan the 
possible values of {rii). From an experimentalist point of 
view, we can apply a varying external potential and mea¬ 
sure the possible values of {ui). We will restrict ourself 
to situations in which fp is invertible, i.e., we can write 
Si = this precludes the possibility of phase 

transitions, for which, for a given Si, there could be more 
than one value of (rii). 

The connection with the description of the previous 
section, based on the nonlinear Fokker-Planck equation, 
is given by the fact that (rii) should be equal to rieq in 
X = ai. Then, using (19), we have 


{ui) = e 


- ^-I3v-s, 


(24) 


V. HARD CORE INTERACTION 

Flard core interaction is the prototypical case study. 
Phenomenological arguments to determine the transition 
probabilities are usually presented [llinilllllH]- With¬ 
out the presence of an external held, the transition prob¬ 
ability from i to i -b 1 is proportional to the quantity 
1 — rii+i, that indicates if the target site is free to be 
occupied by the incoming particle: 

=P(l-n,+i). (28) 

Since the transition probability depends on the target 
site, we have 7 = 1 . From Eq. © we obtain Vi+i = 
-/3“Mn(l - n^+l). In the continuous limit, and neglect¬ 
ing fluctuations, 


In the same way that we considered that the transition 
probabilities satisfy the detailed balance relation even 
when the system is out of equilibrium, we can consider 
that the relation between the mean held potential V and 
the particle concentration that we can derive from the 
previous equation also holds out of equilibrium. This 
can be justihed using the local equilibrium assumption. 
For simplicity, we replace (ui) by n. Therefore, from (24) 
we obtain 


V{n) = -j3 ^ (lnn + fp^{n)Y 


(25) 


The most simple situation is zero interaction; in this 
case n = e~^' and fp^{n) = — Inn. A more direct rela¬ 
tion between V and $ can be derived from (24); it also 
clarihes its physical meaning. We have 


_ 


dlnZi 


_ e ‘ rii 

“ ^ 


— (Sirij —/3$(ni) 


Ui—O 




= -y _ - , 

- 1 )' 

00 

— __ \ ' + 

z, ^ n'f 

n'=0 


= -y 


-/3($(n'-|-l)-3>(n')) 

C / I ^ 

n'! 


n'—O 




(26) 


We can interpret this result using the Jarzynski equal¬ 
ity m, that gives a relation between the Helmholtz free 
energy variation AF and the applied work W: 


= {e-^y. 


(27) 


Comparing (26) and (27), we can see that the mean held 


potential H at a given point is equal to the free energy 
change when the work needed to increase the number of 
particles by one, W = ^{rii -b 1) — 4’(ni), is applied to 
that point. 


V{n) = —j5 ^ln(l —n). (29) 


The resulting equilibrium concentration, that is obtained 
replacing (29) in (19), corresponds to Fermi-Dirac statis¬ 
tics: rieq = 


1 ). 


On the other hand, using the interaction energy 
0 if rii = 0,1 


, \ f 0 if rii = 0 

“ I 00 if rii > 2 


(30) 


in (26), we can arrive to the same result (29). This proce¬ 


dure does not appeal to an a priori dehnition of the tran¬ 
sition probabilities with the shape of (28), and makes a 


clear distinction between the mean held potential V and 
the microscopic interaction energy $. 

Several results obtained from diffusion in a lattice in¬ 
dicate that the hard core interaction does not have any 
effect on the collective diffusion coefficient US mils]: 


D-B.C = Dq. 


(31) 


(Of course, it does have an effect in the single particle 
diffusion coefficient [MIISS].) This result is obtained from 
the expression for the diffusion coefficient (16) using (29) 
for the mean held potential and 7 = 1 . 

Therefore, in our description based on the 7 -dependent 
energy landscape, the result (31) corresponds to a con¬ 


stant value of 7 equal to 1. But this is not the only 
possibility. For example, for 7 = 1/2 and the same hard 
core interaction we obtain 


A’hc = 


Do 

1 — n 


(7 = 1 / 2 ). 


(32) 


It is interesting to note that a similar dependence of the 
diffusion coefficient on the concentration has been ob¬ 
tained for hard spheres in a continuous space (see, for 
example, [a cni HU El]). More complex expressions of 
the diffusion coefficient can be obtained if a 7 parameter 
that depends on the concentration is considered. In all 
cases the interaction, and the equilibrium solution, is the 
same. 

Now we consider the situation in which multiple oc¬ 
cupancy is allowed with a maximum number of particles 
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N. Let us call rrii the number of particles in site i to 
distinguish this case from the one of the previous para¬ 
graphs. And let us consider that the average concentra¬ 
tion is given by (toJ = N{ni). In terms of the partition 
function: 


{rrii) = -N 


dlnZ^ 

~W~ 


ainzf 


(33) 


Then, the partition function for the multiple occupancy 
case is 


N 


Zl = (1 ■ 


3-5. )JV ^ y- g-5.m 

' \ m 

rrii—d 


N 


(34) 


Using the definition of the partition function, and com¬ 
paring with the previous equation, we obtain the inter¬ 
action energy for this case: 


4>'(mj) = 


—/3 ^ In 
oo 


AT! 


(jv-m.)! ifO<m, <1V 
it rrii > N 


The mean field potential is 

V{m) = —/3~^ ln(A^ — m). 


(35) 


(36) 


By increasing the value of N we increase the mean num¬ 
ber of particles in each site and reduce fluctuations, a 
procedure that favors the conditions for the validity of 
the Fokker-Planck equation. 

Figurej^shows the equilibrium solution with a constant 
force to the right in a closed system, and Monte Carlo 
simulation results for different values of 7 ( 0 , 1/2 and 1 ); 
as expected, in all cases the same equilibrium solution is 
obtained. 

As mentioned before, the different values of 7 be¬ 
come relevant in non-equilibrium situations. We ana¬ 
lyzed two simple cases. The first one is a constant force 
{U = —Fx) applied to a system with periodic bound¬ 
ary conditions. The non-equilibrium stationary state is 
homogeneous with non-zero current. From Eq. (|14[ ) we 
can obtain the mobility: /r = Ja/{nF). Figure ^ shows 
numerical simulation results of the mobility against con¬ 
centration for different values of 7 . The results coincide 
with the analytical curve (151. The analytical expres¬ 


sions of /r and D, for fermions and bosons (next section) 
and for different values of 7 , are shown in table |T] 

The second non-equilibrium situation that we consid¬ 
ered is zero force with unequal fixed boundary conditions. 
The difference in particle concentration at both ends of 
the system produces a constant current in the station¬ 
ary state. Numerically, we can obtain the density pro¬ 
file and also its space derivative. From the equation for 


the current (14), we have that the diffusion coefficient is 


D = —Ja/ In this way, we can plot the diffusion co¬ 
efficient against concentration for different values of 7 , as 
shown in hgure|^ The numerical results, again, coincide 
with the analytical results (16) (see table |l]). 



FIG. 2. Equilibrium concentration n for hard core interaction 
with U{x) = —Fx, PFa = 0.05 and total concentration 0.5. 
The curve corresponds to the Fermi-Dirac distribution and 
the dots to Monte Carlo simulations with different values of 7 : 
7 = 0 (circles), 7 = 1/2 (squares) and 7=1 (diamonds). The 
same distribution is attained for all values of 7. Simulations 
were made with multiple occupancy of sites: rii = rrii/N. 
Number of samples: 1000; MC steps for each sample: 10^ 
(7 = 0 and 1, A = 100) or 10® (7 = 1/2, N = 500); fixed 
boundary conditions; system size L = 100a. 



FIG. 3. Fermion’s mobility /r against concentration n for 
different values of 7, obtained from a system with homoge¬ 
neous concentration, constant force PFa = 0.05 and periodic 
boundary conditions (see text for more details). The curves 
correspond to Eq. (151. Dots correspond to Monte Carlo sim¬ 
ulations with 7 = 0 (circles), 7 = 1/2 (squares) and 7 = 1 
(diamonds). Number of samples: 1000; MC steps: between 
10® and 10®; N = 100; system size L = 100a. 
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FIG. 4. Fermion’s diffusion coefficient D against concentra¬ 
tion n for different values of 7, obtained from a system with 
unequal fixed conditions at the ends and zero force, in a non¬ 
equilibrium stationary state (constant current^^ see text for 
more details). The curves correspond to Eq. (16 1 . Dots cor¬ 
respond to Monte Carlo simulations with 7 = 0 (circles), 
7 = 1/2 (squares) and 7 = 1 (diamonds). Number of samples: 
1000; MC steps: between 10* and 10®; N = 500; system size 
L — 100a; n(0) = 0.8 and n{L) = 0. 


TABLE I. Mobility and diffusion coefficient for different val¬ 
ues of 7, as derived from equations (151 and ( 161 , for fermions 
(minus sign) and bosons (plus sign). 


7 

f^/{DoP) 

D/Do 

0 

1/(1 ±n) 

l/ilFnf 

1/2 

1 

1/(1 ±n) 

1 

1 ± n 

1 




FIG. 6 . Boson’s mobility fj, against concentration n for dif¬ 
ferent values of 7, obtained from a system with homogeneous 
concentration, constant force PFa = 0.05 and periodic bound¬ 
ary conditions. The curves correspond to Eq. (15l. Dots 
correspond to Monte Carlo simulations with 7 = 0 (circles), 
7 = 1/2 (squares) and 7 = 1 (diamonds). Number of samples: 
5000; MC steps: between 10* and 6 10®; N = 20; system size 
L = 100a. 


VI. BOSONS 

Boson statistics can be obtained in a classical context 
using a statistical potential. The behavior of quantum 
non-interacting particles can be reproduced by classical 
particles with this effective attractive interaction, that is 
given by 


$(ni) =-/3 Mnrii! (37) 


FIG. 5. Equilibrium concentration n for bosons with energy 
U{x) = —F X, pFa — 0.03 and total concentration 0.5. The 
curve corresponds to the Bose-Einstein distribution and the 
dots to Monte Carlo simulations with different values of 7: 
7 = 0 (circles), 7 = 1/2 (squares) and 7=1 (diamonds). In 
order to reduce simulation fluctuations, we define the concen¬ 
tration Ui at site i as rii = nXi/N, where rrii is the number of 
particles in this site. Number of samples: 1000; MC steps for 
each sample: 10^; N = 10; fixed boundary conditions; system 
size L — 100a. 


Its effect is to cancel the Gibbs factor in the partition 
function (22), from which the Bose-Einstein distribution 


is obtained. 

It is not difficult to obtain the corresponding mean field 
potential from (26): 


V{n) = —P ^ln(I-l-n). 


(38) 


The equilibrium Bose-Einstein distribution is recovered 
when this result is replaced in (19): rieq = — 

1). Figure]^ shows this equilibrium solution with a con- 
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FIG. 7. Boson’s diffusion coefficient D against concentra¬ 
tion n for different values of 7 , obtained from a system with 
unequal fixed conditions at the ends and zero force, in a non¬ 
equilibrium stationary state. The curves correspond to Eq. 
( |16| |. Dots correspond to Monte Carlo simulations with 7 = 0 
(circles), 7 = 1/2 (squares) and 7 = 1 (diamonds). Number 
of samples: 5000; MC steps: between 10^ and 10®; N = 100; 
system size L = 100a; n(0) = 3.5 and n{L) = 0. 


slant force to the right. The figure also presents numer¬ 
ical results for different values of 7 , showing that in all 
cases the same solution is attained. 

As in the case of hard core interaction (or fermions), 
the mean field potential determines the equilibrium so¬ 
lution but not the dynamics. The mobility and diffusion 
coefficient are not unequivocally determined by V^ they 
depend also on 7 . As in the previous section, we analyzed 
two non-equilibrium steady state situations from which 
the mobility and the diffusion coefficient against concen¬ 
tration can be obtained. Numerical results coincide with 
the corresponding equations (151 and (16), as shown in 
figures and 

In all cases, for fermions or bosons, in the limit of low 
concentration, both the mobility and the diffusion coef¬ 
ficient coincide with the corresponding non interacting 
values: and Dq respectively. 


VII. GENERALIZATIONS 


The generalization of the nonlinear Fokker-Planck 
equation (18) to higher dimensions is straightforward: 


dn 

Ik 


S/ ■ n + D\/n). 


(39) 


The expressions for /r and D, (15) and (16), remain un¬ 
changed. The relation between V and $, (26), also keeps 
its validity for higher dimensions. Let us note that these 
equations also hold for a space dependent 7 parameter. 

For the generalization to continuous systems we have 
to take into account some considerations. The system is 
divided in cells of size a, small enough to be considered 
point-like and, at the same time, large enough to contain 
many particles (a standard approach in non-equilibrium 
statistical mechanics). The cell size is much smaller than 
a typical concentration wavenumber A. Also, the inter¬ 
action range r should be much smaller than the cell size. 
So, we have r ^ a A. 


Condition a <C A allows us to take the continuous limit 
and obtain the Fokker-Planck equation. And condition 
r a allows us to neglect the interaction between cells, 
since the interaction energy in the cell’s surface is much 
smaller than in the bulk, and to keep the validity of the 
relation ( |^ between V and $. The consequence is that 
the Fokker-Planck equation is local. 


VIII. CONCLUSIONS 


We derived a non-linear Fokker-Planck equation for in¬ 
teracting particles. The derivation is based on an energy 
landscape of wells and barriers on a lattice. The mean 
field potential and the external potential determine the 
energy wells and the equilibrium solution. The barrier’s 
heights depend on the 7 parameter; it determines if the 
transition probability depends on the mean field poten¬ 
tial of the origin site, the target site, or on a mixture of 
both. 

A relation between the mean field potential and the mi¬ 
croscopic interaction energy was deduced. The Jarzynski 
equality can be used to interpret the mean field poten¬ 
tial as the free energy change when the work needed to 
increase the number of particles by one is applied. 

The results are illustrated with the hard core (or 
fermion) and boson interactions. The corresponding 
mean field potentials can be combined with different val¬ 
ues of 7 . The consequence is that, for the same potential 
and equilibrium solution, the dependence of the mobil¬ 
ity and diffusion coefficient on the concentration can have 
large variations determined by the value of 7 . In all cases, 
for small concentration, mobility and diffusion coefficient 
tend to the corresponding non interacting values. 

We considered only different constant values of 7 but, 
in general, it could be a function of the position or of the 
concentration. 
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